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ABSTRACT 



We propose a method to improve lattice operators composed of Wilson fermions 
which allows the removal of all corrections of 0(a), including those proportional 
to the quark mass, leaving only errors of 0(a 2 ). The method exploits the fact 
that chiral symmetry is restored at short distances. By imposing this requirement 
on correlation functions of improved lattice operators at short distances, the 
coefficients which appear in these operators can be determined. The method is 
an extension of the improvement program of the Alpha collaboration, which, 
up to now, has only been applicable in the chiral limit. The extension to quarks 
with non-zero masses is particularly important for applications in heavy quark 
physics. 



1 Introduction 



An important source of errors in present lattice computations is the use of a finite lattice 
spacing, a. A systematic method for reducing discretization errors order by order in a was 
proposed by Symanzik |] and developed in ref. 0. It consists of modifying the action and 
operators by "irrelevant" terms, chosen in such a way that the convergence to the continuum 
limit is accelerated. For Wilson fermions, in the first implementations of this procedure, the 
improvement coefficients were computed in perturbation theory, leaving errors of 0(ag^) in 
physical quantities 0, |J, where g is the bare coupling constant. It turns out, however, 
that, for lattice spacings used in present simulations, one-loop perturbation theory is not 
sufficiently accurate for some quantities |5], even when "tadpole resummed" [[?]]. Recently, 
a method for determining the improvement coefficients beyond perturbation theory has been 
proposed and implemented, thus achieving full 0(a) improvement || |J Q This method 
allows one to obtain the improved action and the improved vector currents, but can only 
be applied to other composite operators, including the axial currents, in the chiral limit. 
We explain the reasons for this limitation in section |3]. It is particularly significant for 
applications involving heavy quarks, where 0(am) errors may be large 0. Important examples 
in the phenomenology of .B-mesons include the calculations of the leptonic decay constant 
fs-i the form factors of semileptonic B decays, the S-parameters and the amplitudes of the 
radiative decay B — > K*j. 

In this paper we suggest a new method, based on the short-distance behaviour of gauge- 
invariant correlation functions, which overcomes this limitation. The basic idea is to exploit 
the restoration of chiral symmetry at short distances in the renormalized theory. This is 
achieved by constructing finite improved operators by matching lattice and continuum cor- 
relators at short distances. The matching is implemented by requiring the vanishing of 
chirality violating form factors in the correlation functions. This procedure is justified by 
the following two observations. Firstly, at short distances, renormalized perturbation theory 
becomes chirally invariant (explicit chiral symmetry breaking effects induced by the regular- 
ization are reabsorbed by imposing the validity of the Ward Identities of chiral symmetry, 
while violations from the non- vanishing of quark masses disappear at short distances). Sec- 
ondly, contributions due to the spontaneous breaking of chiral symmetry, which are absent 
in perturbation theory, vanish in this region. So both effects decrease as we go to short dis- 
tances. In all the cases of interest, including the matrix elements of quark bilinear operators 
considered in this letter, the terms which violate chirality vanish either quadratically with 
the quark mass, i.e. as m 2 |x| 2 , where \x\ is the Euclidean distance, or linearly as mAQ CD |a;| 4 . 
On the other hand the lattice artefacts vanish as ma. It is the different short- distance be- 
haviour of the 0(ma) terms which allows us to eliminate them. Note that in this paper 
we are proposing a method for the determination of the coefficients of the 0(ma) terms. 
For this reason the quark mass is a free parameter, which can be chosen to optimise the 
procedure and does not necessarily coincide with any physical quark mass. For example we 
may envisage using a mass of the order of a few hundred MeV to extract the coefficients, 
and use these coefficients in studies of heavy quark physics. 

The reason why mass corrections, other than the lattice artefacts, vanish either as m 2 \x\ 2 
or as mAQ CD |x| 4 , can be easily understood. The chiral structure of QCD is such that 

1 By full 0(a) improvement we mean that the remaining discretization errors are of 0(a 2 ). 

2 By m we mean some choice of physical quark mass. We do not need to pick a specific definition in this 
paper. 



contributions to gauge invariant correlation functions which arise in perturbation theory must 
contain an even number of mass insertions. For this reason the short distance perturbative 
corrections vanish at least as m 2 \x\ 2 . On the other hand contributions due to spontaneous 
symmetry breaking require the insertion of the mass operator, mqq. By dimensional counting 
this implies that the corrections are of the form mAQ CD |x| 4 . 

Before presenting the technical details of the method, we explain what we mean by "short 
distances" throughout this paper. In the continuum we simply mean that |x|Aqcd «C 1, 
where \x\ is the Euclidean distance. On the lattice, in addition, we require that \x\/a is 
sufficiently large so as to avoid lattice artefacts (a is the lattice spacing). The existence 
of a window, such that a <C \x\ <C Aq^ d , requires that we have a sufficiently small lattice 
spacing. 

To implement a full 0(a) improvement away from the chiral limit, further conditions should 
be satisfied. Indeed, it is necessary to vary the lattice spacing a, distance |x| and quark mass 
m in such a way that the chirality violating terms of 0(m 2 \x\ 2 ) and 0(mAQ CD |a;| 4 ) are small 
enough to isolate the 0{ma) artefacts in the correlation functions. We postpone the detailed 



description of these conditions until subsection where a specific example is discussed. 
For the remainder of the paper we will assume that \x\ and m are chosen so that the terms 
of 0(m 2 \x\ 2 ) are much larger than those of 0(mAQ CD |x| 4 ) (this condition is straightforward 
to satisfy provided the lattice spacing is sufficiently small). 

We wish to stress that the requirement of the existence of the window a <C |x| <C Aq^ d is 
common to all renormalization schemes implemented using computations at a single value 
of (3 (=6/(7o) ■ 111 principle the constraint \x\ <C Aq^ d can be avoided by using a sequence 
of lattices of decreasing a, and correspondingly of decreasing physical volumes. The idea 
is to determine the normalization constants by imposing renormalization conditions on a 
lattice with a very small spacing. On such a fine-grained lattice, very small distances \x\ 
are possible (so the renormalization scale, in momentum space, can be chosen to be very 
large) but the small physical size precludes the calculation of hadronic correlation functions. 
It is therefore necessary to determine the equivalent normalization constants for the larger, 
but coarser, lattice, on which physical matrix elements are computed. This is achieved 
by increasing the lattice spacing a and matching the normalization constants at fixed 
and then by increasing \x\ and matching at fixed a. This procedure is repeated until the 
required value of a is reached. This method has been proposed within the context of the 
Schrodinger functional approach in ref. ||, but can also be applied here, and is illustrated 
with an example in the appendix. 



In ref. JTO] a method of non-perturbative renormalization was proposed, based on the im- 
position of renormalization conditions on Green functions of operators between quark and 
gluon states in a fixed gauge. It is possible to generalize this procedure to implement im- 
provement, exploiting the chiral behaviour of Green functions at large momenta, but this 
adds a number of theoretical complications. In a fixed gauge the basis of higher dimensional 
operators which must be added to the bare lattice ones is larger than with the gauge in- 
variant method proposed here, particularly when computing off-shell Green functions. This 
basis is restricted by BRST-invariance rather than gauge invariance, and includes operators 
which vanish by the equations of motion, but which are not themselves BRST-invariant Jll| . 



A detailed description of the procedure applied to quark bilinear operators will be presented 
in ref. |Tl|, together with an exploratory computation of all the subtraction coefficients. 



The use of chiral Ward identities to impose constraints on the normalization constants of op- 



erators is a well established technique [0 . In particular, for Wilson and tree- level improved 
theories, they provide the normalizations of the vector and axial-vector currents, the ratio 
of the normalizations of the scalar and pseudoscalar densities and the values of the finite 
coefficients describing the mixing of lattice operators of different chiralities (this mixing is 
present due to the breaking of chirality by the Wilson term). Chiral Ward identities also 
provide important constraints for the determination of the operators necessary for full 0(a) 
improvement [^, £| and, where possible, we propose to exploit them to determine some of the 
coefficients. As explained in some detail in section |3|, however, it is not possible to determine 
all the improvement coefficients by the use of the Ward identities only. 

One can also attempt to construct an improved action and operators by plotting the be- 
haviour of the matrix elements of these operators as a function of the lattice spacing, and 
fixing the improvement coefficients by requiring that there are no linear terms in a (including 
possible logarithmic factors). Such a procedure is impracticable, however, at least without 
a considerable amount of additional input based on the techniques described in refs. || 
and in this letter. 

The plan of the remainder of this letter is as follows. In the next section we briefly review 
the construction of the improved action and the determination of the mass dependence of 
the bare coupling constant in the unquenched theory. In section |3| we explain why the 
chiral Ward identities cannot be used to determine all the coefficients needed to construct 
improved lattice composite operators. The description of our method, and its application 
in the construction of improved lattice quark bilinear operators for quarks with non-zero 
masses is presented in section [|. Finally section |5| contains our conclusions and a discussion 
of future prospects. 



2 Improvement of the action 

The Wilson fermion action (Swuson) * s improved by the addition of the SW operator |3j, 

S = Sgimge + ^Wilson + a J (l^X O SW (1) 

where S ga uge is the gauge action and 

Osw = -^c sw ^qcr^F^q , (2) 

with <j^ u = (i/2)[ r y tl , 7„]. Here and in the following we use continuum notation to refer to 
lattice quantities. A method for the non-perturbative determination of csw was presented 
and exploited in refs. 0, || . 

In order that physical quantities should approach their continuum values with deviations of 
only 0(a 2 ), we need, in addition to choosing the correct value of csw, to adjust the bare 
gauge coupling go and the bare mass mo in a way which depends on the renormalized quark 
mass m ||. To see this, consider the renormalized coupling constant g(fi 2 ) defined, for 
example, as the three gluon vertex (in some fixed gauge) at large virtualities characterized 
by the scale /i 2 . The corresponding renormalization constant Z g , defined by 



g{n 2 ) = Z g g (a) , 



(3) 



contains terms of 0(ma) from diagrams which include quark loops, which must therefore 
be removed for improvement of unquenched lattice QCD. Following ref. || we absorb this 
linear term into a redefinition of the bare coupling constant go, 

9o =9o0- + b g ma) , (4) 
so that Z g has no remaining linear term in ma. 

Thus, for each value of the bare quark mass, we need to determine go []. For simplicity 
we take the quarks to be degenerate, so that there is only one mass parameter mo; the 
generalization of the following discussion to the non-degenerate case is straightforward. We 
can imagine, at least in principle, proceeding as follows: 

• Choose a value of go, compute a short-distance physical quantity Q at several values 
of the quark mass at this value of go and extrapolate the results to the chiral limit. An 
example of a suitable quantity is the derivative of the potential between a static quark 
and antiquark (V'(r)) at a separation of a fixed number of lattice spacings (r = na), 
satisfying tAqcd <C 1 and rm <C 1. We denote the result by Q(n,go,0), where the 
final argument represents the fact that this value corresponds to m = 0. 

• The quantity Q must be such that its linear dependence on the mass is entirely due to 
lattice artefacts, and that these 0(ma) effects can be removed by fixing csw and b g only. 
For physical quantities, such as decay constants and form factors, which are related 
to matrix elements of composite operators, there are further mass dependent terms 
of 0(ma) in their renormalization constants. It is therefore not possible in general to 
use these quantities to determine b g (with the exception of the matrix elements of the 
vector current, see subsection fO| below) . 

• We now choose a value of mo and recompute Q at different values of the bare coupling 
until we obtain the same result, i.e. until we find a coupling g such that Q(n, go, m) = 
Q{n, go, 0). In this way we obtain g for each mo- As discussed in the introduction the 
result is independent of n, up to quadratic corrections of order (mr) 2 . 

• Finally, we need to determine the values of the lattice spacing corresponding to go{ a ) 
and of the bare quark mass (m ) of the physical quark. This can be done by simu- 
lating with a series of values of {m ,^o( m o)} and requiring that the computed results 
for two physical quantities, for example the masses of two hadrons, agree with their 
experimental values. 



The technique described above needs considerable computing resources, since it requires 
new simulations at a number of different values of (3 in order to map out the behaviour of 
go with mo- At this stage we are not very concerned about this, since it will be some time 
before an accurate determination of b g is actually required. It will only be needed when 
precise unquenched computations become possible (b g can be taken to be zero in quenched 
calculations). Even then, however, the effects of improvement (beyond e.g. tree- level) will 
be most significant for heavy quarks, for which the effects of quark loops are expected to 
be strongly suppressed. We also note that in first order perturbation theory b g is small, 
b g ~ 0.012 Nf g>Q, where Nf is the number of quark flavours |13| 0. Thus, for the purposes of 

3 See ref. [|| for a method based on the use of Schrodinger Functionals. 

4 The one-loop results for other improvement coefficients can be found in ref. 1 14 



this letter, we simply note that a non-perturbative computation of b g is possible in principle 
and postpone any detailed consideration of the most practical way of determining it. We 
focus instead on the more immediate problem of determining improved composite operators, 
thus, in particular, removing the 0(ma) effects due to heavy valence- quarks. 

For the remainder of this letter we assume that csw is known, and that, as the bare quark 
mass m is varied, the bare coupling constant has been readjusted to the corresponding value 
go(m ). 

3 Chiral Ward identities 

Before describing our method for the construction of improved lattice composite operators, 
we explain why it is not possible to use chiral Ward identities to obtain all the necessary 
coefficients. As an example, consider a lattice quark bilinear operator qTq, where F is one 
of the Dirac matrices, and we have suppressed the flavour indices. The improved lattice 
operator takes the form 

O r = Z r (qTq + ac r 4 ,r) , (5) 

where the O^r are possible dimension 4 operators with the same quantum numbers as Or 
(which will be exhibited explicitly in the following section; there are no such operators for 
the scalar or pseudoscalar densities). The operators on the right-hand side of eq.(|5|) are 
the bare lattice operators. O is constructed so that its matrix elements between physical 
states have discretization errors only of 0(a 2 ). Following ref. H, we parametrize the mass 
dependence in the overall renormalization constants by 

Z v (m) = Z r (0)(l + b r ma) . (6) 

The point which we wish to stress here is that the Ward identities do not determine the b^s. 
In the remainder of this section we briefly explain the reasons for this limitation and in the 
following one we propose a new method for the determination of the &r's- 

Consider the following continuum Ward identity 

(^Jj±x[2mP{x) - d^x)] 0(y)-50(y)} 1 (z 1 ) 2 (z 2 ) ■ ■ ^ = , (7) 

where P and are a pseudoscalar density and axial current respectively, and O, Oi, 2 etc. 
are composite operators, which, for the purposes of this letter, we can assume to be gauge 
invariant. We are free to make any convenient choice of the operators {Oi,02, ■ • •}. For 
simplicity of presentation we suppress the flavour indices on all the operators. The integral 
is over a region R which includes the point y, but excludes the points Z{. SO is the axial 
variation of O. The improved lattice version of the operator 2mP(x) — d^A^x) is obtained 
(up to an overall normalization) by requiring it to satisfy Ward identities of the form (up to 
terms of 0(a 2 )) 

([2mP(x)-d fl A fl (x)]0 1 (z 1 )0 2 (z 2 )---) =0, (8) 

where x is separated from the other points Zi. It might therefore be hoped that by subse- 
quently imposing the identity in eq.(^) it would be possible to determine the improved forms 
of O and 50 (up to an overall normalization factor). This is not so, however, since in order 
to satisfy this identity up to terms of 0(a 2 ), using lattice operators, we must add irrelevant 



operators to P(x) and 0(y) which "vanish" by the equations of motion, and which therefore 
do not contribute to physical matrix elements. For example, the improved form of Or m 
off-shell matrix elements takes the form 

O r = Z T (qTq + ac' r q T($ + m ) + (-^ + m )T q + ac r 4 ,r) , (9) 

where + tjiq is the lattice version of the fermion matrix appearing in eq. ([!]). The terms 
proportional to the new coefficients c' r contribute to the identity in eq. (|7]) in the region where 
P(x) comes into contact with 0(y) (]. Thus, because of these "contact terms" , there are extra 
coefficients, the c r 's, to determine, and in fact it is not possible to obtain the c r 's and the 
o r 's separately using the Ward identities (and so it is not possible to obtain the o r 's which 
are needed for the evaluation of physical matrix elements). 

We now justify this claim in more detail. Using the expressions ([5]) for the operators P, A^, O 
and SO in eq.(|7|), we find that the term containing the rotated lattice operator 50 is 



ZaZq 



2pa(c' P + c' ) 



{50{y)O l {z 1 )0 2 {z 2 )---) , (10) 



ZaZo 

where p is computed directly from the lattice correlation functions in eq.(|8|): 

_ (d,A,(x)Q 1 (z 1 )Q 2 (z 2 )...) 
9 (P(x)O l (z 1 )0 2 (z 2 )---} • 1 ' 

In eq.([TT|), x is separated from the {zi}, and p is independent of the choice of the operators 
{Oj} and of all the coordinates. Since the term proportional to p in eq.(|T0|) is multiplied 
by a factor of a, it is sufficient to take any convenient, not necessarily improved, lattice 
pseudoscalar density P or axial current in eq.(|II]) f\. From eq. ( |T0"D we see that the Or's 
and the c r 's appear in the combination: 

Z A t)Zo(0) { 1 + {b5 ° ~ b °~ hA) ma } " 2pa{C ' P + C ' o) ' (12) 

where the fer's are defined in eq.(^|). For each bilinear O, the 0(a) part of the Ward Identity 
in eq.([7D allows one to obtain the combination of 0(a) that appear in eq.(|T^). Thus one only 
obtains a linear combination of the &r's and the c r 's, and even by considering all possible 
bilinears it is not possible to determine any of the &r 's separately. 

Despite this shortcoming, chiral Ward identities should be part of any improvement program; 
they provide a determination of the cr's as well as many checks of the numerical results for 
the remaining improvement coefficients. 



4 Improvement of bilinear quark operators 

In this section we propose a new method for the determination of the improvement coef- 
ficients in general, and the 6r's in particular. We illustrate the method, by discussing the 
construction of improved quark bilinear operators. We consider only flavour non-singlet 
bilinears, but, unless specifically required for the discussion, we drop the flavour indices. 

5 Note that the integral of d^A^ in eq.(Q) gives a surface integral, and so produces no contact terms. 

6 In this letter we have not needed to define an improved quark mass. However, comparing equations 
(||) and (pj]), leads us to write p — tjiZp/Za- By computing p using improved operators and P, whose 
construction is explained in section |J, the corresponding mass m is indeed a renormalized improved quantity. 



4.1 Pseudoscalar and scalar densities 



The renormalized, improved pseudoscalar density has the general form 

P{x) = Z P P(x) , (13) 

where the renormalization constant Zp = Zp(m) = Zp(0)(l + bp ma + • ■ ■) is, in general, a 
function of the quark mass and P(x) = q{x) n f^,q{x) is the bare lattice pseudoscalar density. 
We start by constructing the two-point correlation function Gpp(x) = (P(x)P(O)) and by 
studying it in the limit \x\ — ► 0. By u \x\ — > 0" we mean "short distances" in the sense 
explained in the Introduction. In this region, we expect on the one hand to be able to use 
the operator product expansion, keeping only the leading, most singular term, and yet on 
the other hand to neglect contributions from contact terms. 

At any given value of m, provided ma < 1, we determine Zp by imposing, at short distances, 
the following renormalization condition 

Gpp(x) = Z P (m)(P{x)P(0)) 

= (P(x)P(0)) cont as |ar| — ► 0. (14) 

Gpp(x) ee (P(x)P(O)) is the bare lattice pseudoscalar-pseudoscalar correlation function, 
computed by numerical simulation at the quark mass m. In eq.flTJl) (P(x)P(0)) cont is the 



pseudoscalar-pseudoscalar correlation function computed in perturbation theory in some con- 
tinuum renormalization scheme, e.g. in the MS or the RI schemes [|l0j at a renormalization 
scale /i. A convenient choice is fi ~ l/\x\. At short distances (P(x)P(0)) cont is independent 
of the quark mass, up to subasymptotic corrections vanishing as m 2 \x\ 2 . This shows that we 
can determine the 0(am) correction to Zp(m) using only gauge invariant quantities. 

An alternative method to determine Zp is the following. We define, at short distances, 

d r-i , u ^2 Gpp(x)\ m= o 

R P = (1 + bpma) = — — — — , (15) 

Gp P {x)\ m 

where Gpp(x) | m= o is the bare lattice correlator Gpp extrapolated to the chiral limit. Thus 
a computation of the ratio Rp provides us with a determination of bp. In addition, at short 
distances Rp is independent of x, which provides us with a set of consistency conditions. 

The leading physical chirality- violating correction to Rp(x) is of the form (up to logarithmic 
corrections) m 2 |x| 2 . In order to be able to extract bp we require m 2 |x| 2 <C ma, which at 
fixed a imposes constraints on the allowed values of m and \x\. In practice, for values of 
the lattice spacing currently used in simulations (a -1 = 0(2-4 GeV)), these conditions can 
be implemented without difficulty Q In addition, one can estimate the 0(m 2 \x\ 2 ) terms by 
studying the short- distance dependence of Rp(x) on x at fixed m, and try in this way to 
improve the precision in the determination of bp. These considerations also apply to all the 
bilinears considered below. 

The above is a good example illustrating the necessity of knowing b g (or equivalently, of 
knowing g (a) for each value of the quark mass). The numerator and denominator on the 
right-hand side of eq.flToD have to be evaluated at the same values of the coordinates x (in 

7 At short distances, the terms of 0(mAq CD |a;| 4 ) are naturally negligible compared to those of 0(m 2 \x\ 2 ) 
provided that the mass m is not much smaller than Aqcd- 



physical units). As explained in section |2], in order to satisfy this requirement (up to 0(a 2 ) 
corrections), the bare coupling constant has to be varied as a function of the mass as in 
eq.(§. 

Having determined the mass dependence of Zp(m), Zp(0) can be obtained for example, by 
imposing renormalization conditions on the pseudoscalar density in the chiral limit using the 
Schrodinger functional techniques || or by determining Zp(m) on quark states, as suggested 



in ref. [fLO ], and then by extrapolating the results for Zp(m) to the chiral limit. In the latter 
case we can ignore completely the extra terms coming from operators which vanish by the 
equation of motions, or from non-gauge invariant operators, since the systematic error due 



to these terms disappears as m — ► [11 



The determination of the renormalization constant of the scalar density, Zg, proceeds in 
exactly the same way. The ratio of the renormalization constants Zs(0)/Zp(0), which is a 



finite quantity, can be determined in the standard way using chiral Ward Identities [12 . 



4.2 The Vector Currents 

The improved vector current has the general form 

%{x) = ZyV^x) , (16) 

where 

Vp(x) = g(x)7 /U g(a;) + ac v ^ i9 » [Q( x ) a ^q(x)] , (17) 

i/ 

and the operators on the right-hand side of flTT| ) are the bare lattice ones. In this case, 
improvement requires the determination of two constants (Zy and Cy), one more than for 
the pseudoscalar and scalar densities. The necessity of the term proportional to cy for 
improvement has been pointed out and discussed in refs. ||, [15]. For any choice of the quark 
mass, Zy = Zy(m) = Zy (0)(1 + byam) can be immediately fixed by requiring that the 
forward matrix element of the current is correctly normalized, i.e. that 

(p\V,\p) = Z v (m)(p\V,\p) = 2 Pfl . (18) 

This does not require the knowledge of cy since the forward matrix element of d v \g(x)a liV q(x)] 
is zero. We then determine cy by studying the short- distance behaviour of the correlation 
function 

G vv (x) = (0\V„(x)V v (0)\0) . (19) 
cy is to be chosen to remove the residual dependence of Gyy(x) on ma as \x\ — > 0. 

It is also possible to exploit the chiral Ward identities (in the massless limit) to determine 
cy 0. For example, having determined Zy from the normalization of the charge operator, 
and anticipating that we can obtain the improved axial current as described in subsection [O 
below, the Ward identity (see eq. (0) ) 

fabc f A a (t 2 , x) - %(ti, x) ) A\{y) V° source (0) )=iN f ( Vf(y) V* smlrce (0) ) (20) 



This is the method currently being used by the authors of ref. . 



provides an equation for cy. In eq.(pO|) Vf source (0) is any vector current (not necessarily 



improved), we have assumed that < ti < y < t 2 , and a sum over the spatial index i is 
implied. It has been convenient to exhibit the flavour indices {a, b,c, ■ ■ •} explicitly By 
varying the parameters in eq. (p0|) (ti t 2,y and the form of Vf source (0) ) we obtain a series of 
consistency conditions. 

We stress that when using Ward identities, such as that in eq. (p0[) , we do not require "short 



distances". Indeed for \y\ — > in eq.(|20D, the contribution from the term containing cy 
vanishes (i.e. is of 0(a 2 )), and hence cannot be determined. To see this, note that this term 
is proportional to 

ac v ( q{y)o^T c q{y) V* source (0) ) , (21) 

which vanishes at short distances by chirality (on the lattice the correlation function in 
eq . (pT|) is of 0(a) because of the Wilson term, leading to a contribution of 0(a 2 )). Thus, in 



order to use the Ward identity to determine cy we need to work at distances y > 1/Aq CD , 
where the non-perturbative chiral symmetry breaking terms make the correlation function 
in eq.([H]) non- vanishing. 

If the Ward identities are used to determine cy, then rather than using the independence of 
Gyy(x) of the quark mass to determine cy, we can use it instead to determine go(rn) (b g ). 
Thus, we see that there are many possibilities to (over-) determine the set of parameters 
{b g ,Zy(0),b v ,c v }. 



4.3 The Axial Currents 

The improved form of the axial current is 

A^x) =Z A A lx {x) (22) 

where 

A?(x) = g(x)7 A1 7 5 g(x) + ac A d^ [q(x)^ 5 q(x)] , (23) 



and the operators on the right-hand side of (|23| ) are bare lattice ones. As in the case of 
the vector current, there are two constants (Z A and c A ) to determine, but now there is no 
normalization condition which would allow us to determine Z A without knowledge of c A . 
We start therefore with the determination of c A . Following || we exploit the fact that 
n(x) is proportional to P{x), and compute 

Ra ^ c ^~ (p(x)p(o)) • (24) 

c A is determined, together with cswi by requiring that R A is independent of x (up to terms 
of 0(a 2 ) which we are neglecting). R A is also independent of the form of the pseudoscalar 
operator at the origin (e.g. -P(O) could be replaced by some extended operator in both the 
numerator and denominator of eq. (|2l|) ; see also the discussion around eqs.(|TT|) and (p0|) ). 
This provides further constraints which can be exploited in the simultaneous determination 
of c A and c sw H- 



9 We use a notation in which, for example, P a (x) = -ip(x)j 5 T a ip(x), with the generators of the flavour 
c 

j-abc are £ ne correS p 0nc img structure constants, [t^t 6 ] = if 



symmetry normalised to Tr(r a r b ) = \8 ab . We present the results for an SU(Nf) flavour symmetry, and 



The determination of Za can now be achieved in a number of ways. For example, one can 
start by calculating Za in the chiral limit by imposing a chiral Ward identity ||, using the 
improved operators V and A. We will see below that in this case we do not need to know 
the values of Zy and Cy. For example, consider the identity 

j ahc YJ, P e {z) ( A a (h, x) - A a (t u x) ) 3*(y) P d (0) )=iN f ^2( P e (z) V Q %y) P d (0) ) , 

x,y y 

(25) 

where < ti < y < t 2 < z$ and again it has been convenient to exhibit the flavour indices 
{a, b, c, ■ ■ •} explicitly Any convenient form for the lattice pseudoscalar densities, P, can be 
used in eq. (|ZoT) ; indeed the pseudoscalar densities could also be replaced by other operators. 
Both sides of eq.(p5|) can be simplified. Current conservation implies that the left-hand side 
is independent of tx,2 and yo (provided that the above time-ordering is preserved) and hence 
that the two terms are equal (note that as we are summing over x and y, we can, for example, 
change variables x — > y in the second term). Combining this with the normalization of the 
charge operator on the right-hand side we obtain 

N 2 

f abc f cde Yj pe W 450*) #o(y) P d (v) > = — ^ ( pc W pc (°) ) > (26) 

where < y < x < z Q . For each set of values of {x ,y ,z } satisfying this ordering, and 
for each choice of z, eq. (|26|) is an equation for the one unknown Za{0) (ca having been 
determined as described above). Finally, the values of Za{to) = Z^(0)(1 + bAfna) away from 
the chiral limit can be obtained from the condition that at short distances the correlation 
function (A fl (x)A l ,(0)) is independent of the quark mass, 

(^(i)i(0))| m = {^(x)^(0))U =0 . (27) 

For any value of the quark mass, eq.([27]) represents an equation for 6^. With this method, 
no perturbative calculations, either on the lattice or in the continuum are necessary. 

Alternatively one can calculate the correlation function ( A^(x)A u (0) ) at short distances 
in continuum perturbation theory (here A^ is the continuum axial current and, at short 
distances, the correlation function is independent of the mass) and determine Z^(m) by 
imposing that ( A^(x)A u (0) ) is equal to this perturbative result. Finally, as done for the 
pseudoscalar density, we can compute Za on quark states, for different masses, and obtain 
Za{0) by extrapolating the result to the chiral limit. This gives the correct result provided 
the renormalization scale is large enough and the renormalization conditions are compatible 
with the Ward identities for the axial current |TI5 . 



4.4 The Tensor bilinears 

Finally, we briefly sketch the corresponding analysis for the tensor operator. The improved 
tensor operator has the form 

% v {x) = ZtT^x) (28) 

where 

Tfa,(x) = iq{x)a, xlJ q{x) + ac T (d^V„ - d v V^) , (29) 



and the operators on the right-hand side of fl29|) are bare lattice ones. Since we are neglecting 
terms of 0(a 2 ), we do not need to include the 0(a) terms proportional to cy in the expression 
for V p given in eq. flTT|) . 

Again we have two constants (Zt and ct) to determine and, as in the case of the axial 
current, we start with the determination of ct- This can be achieved by noting that the 
correlation function ( T pv (x) V p (0) ) is zero at short distances. This is because the tensor and 
vector currents have different chiralities, and chirality is a good symmetry at short distances. 
The vanishing of this correlation function represents an equation for ct- 

The coefficient ct can also be determined using Ward Identities in the chiral limit. For 
example, consider the identity 

E E W( ( f ) - 45(*1' f ) ) T ksource(0) > = 2d* t (T%( V ) T^ source (0) ) , 

X \iV 

(30) 

in which no indices are summed unless explicitly indicated, and we have assumed < t\ < 
Do < t<i- This identity is valid for any form of the tensor density at the origin; the simplest 
case is to assume that T^ source = iqa^q. The set of equations fl3U| ) for different choices of 
p, cr, y and T^ VtSOUrce , allow us to determine ct as well as providing many consistency checks. 
Note that Zt appears as an overall factor on both sides of the equation, and so cannot be 
obtained from eq. (p0|) . 

The overall normalization constant Zt can then be determined, for example, by computing 
the correlation function (T^^x) T pa (0) ) at short distances in perturbation theory (in some 
renormalization scheme) and on the lattice and imposing that they are equal. Zj<(m) is now 
the only unknown in this equation and can therefore be obtained. 

Zt can also be obtained in the chiral limit by imposing a normalization condition on the 
matrix element of the tensor operator between quark states. The mass dependence of 
Z^(m) is then determined in the standard way, by requiring that the correlation function 
(T pu (x) T pa (0) ) is independent of the quark mass at short distances. 

5 Conclusions 

In this letter we have proposed a method for the construction of improved lattice composite 
operators, containing quarks with non-zero masses. The remaining discretization errors in 
matrix elements of these operators are of 0(a 2 ). The method builds on the improvement 
program of the Alpha collaboration and extends it away from the chiral limit. This is 
of particular importance in the study of the decays of hadrons containing a heavy quark. 
A detailed numerical study is now necessary to determine the precision with which the 
coefficients can be obtained; the results will be presented in a future paper. 

Our proposed method is based on the observation that chiral symmmetry is restored at short 
distances. By imposing this requirement on correlation functions of composite operators, we 
can determine the coefficients which are needed to construct the improved operators. For 
purposes of illustration, in this letter we have studied the application of our method to 
quark bilinear operators; the approach can, however, also be generalized to deal with more 
complicated composite operators, such as four-quark operators relevant for studies of the 
effective weak Hamiltonian. 



We are also carrying out a pilot study of the determination of the normalization coefficients 
by imposing normalization conditions between quark and gluon states in a fixed gauge. Again 
the key ingredient in the method is the restoration of chiral symmetry at large momenta or 
short-distances. The basis of higher dimensional operators which must be considered (and 
hence the number of coefficients) is larger in this case, and the method will be presented in 



detail in ref. [pi], together with the results of the preliminary numerical study. 
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Appendix 

In this appendix we sketch the procedure which allows, in principle, the evaluation of Zp 
away from the chiral limit, in a way which obviates the need for the condition \x\ <C Aqq D . 
Following ref. ||, we start with a fine lattice with a spacing which is so small that one 
can perturbatively match the lattice results onto a continuum renormalization scheme with 
negligible errors from higher order corrections, and, after successive iterations, we end with 
a course lattice which is sufficiently large to allow the computation of hadronic correlation 
functions. For simplicity we present the discussion in the quenched case only, which obviates 
the need to determine b g . The generalization to the unquenched case introduces further 
technical complications, but the basic strategy is the same. The steps are the following 0: 

i) On a lattice with a very small spacing, ai, and bare coupling gi, we use our method to 
determine Zp(gi, /i, mi) for a given value of the quark mass m\, and for a renormaliza- 
tion scale Aqcd <C [i <C %~ . This condition is readily satisfied since we are working 
on a very fine lattice. We do this using two-point functions at distances ~ 1///. 

ii) We evaluate Rp(x 2 , 9i, m i) as defined in eq. (|l"5|) and a second physical quantity 
Q(x 2 , gi, mi) for which the 0(a) effects can be removed by adjusting csw only (see 
section ||). We choose \x 2 \ to be larger than \xi\, e.g. a simple choice is |a?2 1 = 2|xi|. 

iii) The next step is to obtain the same physical quark mass on both lattices. With a 
coupling constant g 2 , corresponding to a lattice spacing a 2 larger than ai, e.g. a 2 = 2a±, 
we adjust the mass of the quark until Q(x 2 , gi,mi) = Q(x 2 , g 2 ,m 2 ) (x 2 is the same in 
physical units on both lattices). The calibration of the lattice spacing is performed, as 
usual, in the chiral limit, using, for example, Q(x 2 , gi,0) = Q(x 2 , g 2 , 0). 

10 In this appendix, we exhibit explicitly the relevant arguments of Zp, Gpp, Rp etc. where needed. 
Some of these arguments were implicit in the earlier sections. 



iv) We find the new value of bp by imposing the condition Rp(x2, #2, ^2) = Rp(%2, gi, Tax). 

v) At this point we only require Zp in the chiral limit. This is obtained by imposing the 
condition Z P (g 2 , yU, 0)G PP (x, g 2 , 0) = Z P (g 1 , n, 0)G PP (x, g 1 , 0). 

By iterating the steps ii) - v) we can reach a large lattice (in physical units). 

For the other bilinears, one first determines the subtraction constants, cy, ca and op in 
the chiral limit, and then use the above procedure to obtain the overall renormalization 
constants. 
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